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Abstract. We study models in which a light scalar dark energy particle couples to 
the gauge fields of the electroweak force, the photon, Z, and bosons. Our analysis 
applies to a large class of interacting dark energy models, including those in which 
the dark energy mass can be adjusted to evade fifth-force bounds by the so-called 
"chameleon" mechanism. We conclude that — with the usual choice of Higgs sector — 
electroweak precision observables are screened from the indirect effects of dark energy, 
making such corrections effectively unobservable at present-day colliders, and limiting 
the dark energy discovery potential of any future International Linear Collider. We 
show that a similar screening effect applies to processes mediated by flavour-changing 
neutral currents, which can be traced to the Glashow-Iliopoulos-Maiani mechanism. 
However, Higgs boson production at the Large Hadron Collider via weak boson fusion 
may receive observable corrections. 
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1. Introduction 

The emergence of cosmology as a data-driven science in the late 1990s enabled 
our theories of the universe to be promoted from mostly speculation to meaningful 
quantitative investigation. Although many components of what now forms the standard 
"concordance" ACDM cosmology had been proposed prior to the quantitative revolution 
and were found to be consistent with experiment, among the more surprising revelations 
was the emergence of a new scale at around 10~^ eV, associated with an apparent 
acceleration of the cosmological expansion. The properties of Nature at this scale have 
been accessible since the earliest days of particle physics, and our models of microscopic 
processes at these energies are now very well tested. It was therefore surprising to 
discover that this hitherto mundane scale was to be associated with an exotic species of 
matter with energy density A ~ (lO'^eV)"^ and equation of state p ~ —p. 

The microphysics associated with this energy density remains unknown. The 
most parsimonious interpretation of the data requires only Einstein's "cosmological 
term," which we now know to be degenerate with the aggregate effect of quantum 
vacuum fluctations. At a mass scale M these contribute a cosmological effect of order 
M'^. Consequently, if we take our well-tested quantum-mechanical theories of physics 
seriously even at comparatively modest scales (up to Mew ~ 100 GeV - 1 TeV) then 
without a remarkable cancellation we encounter a serious disagreement with the data. 
An alternative interpretation is to suppose that unknown physics renders the quantum 
zero-point energy negligible or unobservable. If this is the case, it is possible that our 
present phase of acceleration is driven by the potential energy associated with some 
scalar field. This field would have to be very light on large scales in our present vacuum, 
with mass of order Hq ~ 10"'^^ eV, but it might evade the stringent bounds associated 
with long-range forces mediated by light bosons if its mass could be adjusted to be large 
in regions of high average density. Theories of this type were proposed by Khoury & 
Weltman |T1 [2], who called such fields "chameleonic" in view of their ability to vary 
their properties depending on the environment I^] 

The chameleon property means that models involving these fields can give rise 
to successful acceleration at late times [5], while remaining consistent with known 
constraints on long-range physics. Such models are attractive for another reason, 
because the requirement that the field can respond to local variations in the density 
of bulk matter means that couplings to Standard Model states are mandatory. 
Chameleonic fields are therefore constrained by precision measurements of the early 
universe — in particular, observations of Big Bang Nucleosynthesis ('BBN') and the 
redshift of recombination [5], [6] . As the universe cools the background dark energy field 
remains fixed in the minimum of its potential, whose location slowly drifts. The result 
is a variation in the mass of any particle to which dark energy is coupled. However, 
acceptable models are constructed in such a way that only small changes in particle 
mass can be expected, and therefore the constraints from observations such as BBN are 

f For earlier work, see Refs. Ollj. 
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rather weak. Interesting bounds have also been obtained from a variety of astrophysical 
and terrestrial processes [7l[8l[9l[T0l[IIl[T2l[l3l|lll[T5]. 

These couplings also imply the existence of an interesting collider phenomenology. 
With the aim of complementing the cosmological and astrophysical tests, our purpose 
in this paper is to take the first steps towards understanding the implications of dark 
energy corrections for Standard Model processes which can be observed at present and 
future particle colliders. A related study has been performed by Kleban & Rabadan 

m- 

What form would these corrections take? We expect that the dark energy scalar 
is not charged under any of the usual gauge quantum numbers associated with the 
Standard Model. Its couplings to Standard Model states are therefore unrestricted by 
considerations of gauge invariance. Nevertheless, because bulk mass in the macroscopic 
world is dominated by hadrons it seems unavoidable for a chameleonic scalar to couple 
to those degrees of freedom charged under QCD, namely the quarks and gluons. 
Unfortunately, hadron interactions in QCD are non-perturbative in nature and are 
difficult to study. It is less obvious that the dark energy is obliged to couple to degrees 
of freedom charged under the electroweak SU{2) x U{1) gauge symmetry, but if it does 
then one might imagine that such interactions would offer a more tractable probe of 
the theory than the complicated colour physics of QCD. Our purpose in this work is to 
study the comparatively clean experimental signatures which arise at low energy from 
the existence of couplings between dark energy and Standard Model states which carry 
electroweak quantum numbers. 

Interactions between a scalar dark energy species and the electroweak sector need 
not be harmless. For example, variation in the dark energy vacuum expectation value 
could lead to a shifting fine-structure constant or loss of conservation of electric charge 
p!7t [18] . From the perspective of collider phenomenology, there is another serious 
difficulty: fundamental scalar fields are well-known to depend sensitively on the details of 
physics in the ultra-violet. If Standard Model particles can radiate into light chameleon 
states while participating in some measurable process, then we must allow for the 
possibility of significant corrections to observable Standard Model reactions. Indeed, it 
is a serious question whether any dark energy model of this type can be compatible with 
existing data. It is also important to understand whether we should expect dramatic 
signals at impending high-precision experiments such as the Large Hadron Collider 
(LHC) at the European Organization for Nuclear Research (CERN) or at a proposed 
future International Linear Collider. 

In this paper, we study the effect of such radiative corrections. Our results 
apply to models of chameleon dark energy, and also to alternatives such as coupled 
quintessence, or any beyond-the-Standard-Model scalar species which is light in the 
laboratory environment. Similar issues have been addressed previously by Einhorn 
& Wudka [19], who determined the criteria for heavy scalar particles to be screened. 
However, our results are not contained in their analysis because the scalar particles 
which can cause successful cosmological acceleration must ordinarily be very light in 
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the laboratory environment compared to the electroweak scale, with masses of order 
< 10~^ eV or lighter. It is the effect of highly suppressed couplings in the laboratory 
which allows such particles to have evaded detection, rather than the significant energy 
cost of producing them in collisions. 

In §2. II we give a brief summary of ultra-violet effects in scalar field theories, before 
going on to review the formalism used to study corrections to electroweak precision 
observables In §3.11 we study corrections to the width of the Z boson (a tree- level 
effect), and show that it leads to a weak constraint. In §3.21 we identify a class of loop 
effects which lead to stronger constraints, the so-called "oblique" corrections. The key 
quantities we require to compute them are the vacuum polarizations of the W^, Z and 
7 bosons, which are obtained in §H In §4.21 we interpret these vacuum polarizations in 
terms of an effective Lagrangian which makes their physical content transparent. In S|5] 
we discuss our findings and indicate how our results could be extended to a larger zoology 
of processes, including so-called flavour-changing neutral currents. In particular, in §5.11 
we discuss the conditions under which the largest corrections are "screened," meaning 
that they do not enter in any measurable relationship between observables. It is only 
when screening occurs that the model is automatically compatible with the simplest 
predictions of the Standard Model. In §5.21 we determine the constraints which can 
be obtained from data obtained by present-day colliders, and discuss the role of future 
hadron-hadron or e~^e~ colliders. Finally, in ^we state our conclusions. Some technical 
details are collected in two Appendices. 

We choose units throughout such that h = c = 1. Our metric convention is 
(—,+,+,+), so that on-shell particles have negative invariant momenta. Spacetime 
indices are denoted by lower-case Latin indices {a, 6, c, . . .}, and we label the species of 
vector bosons by upper-case indices {A, B,C, . . .}. 

2. Electroweakly interacting dark energy 
2.1. Ultra-violet effects 

The problem of sensitivity to ultra-violet effects is universal in any theory of scalar fields. 
While it is an obstacle for model-building, UV sensitivity can be exploited as a tool to 
probe the theory at energies much higher than those which can physically be realized 
in particle accelerators. An important example of this occurs in the Higgs sector of the 
Standard Model, which has many parallels with the case of interacting dark energy. For 
this reason we digress to give a brief discussion of the Higgs case, before returning to 
dark energy in §2.21 

All particles which gain their mass via the Higgs mechanism are entitled to radiate 
into Higgs states, and in consequence it was pointed out long ago by Veltman that 
electroweak quantities can receive large Higgs contributions, up to some scale above 
which radiation is suppressed. This scale is presumably determined by a more complete 
theory of microscopic interactions, of which the Standard Model is an effective low 



Collider constraints on interactions of dark energy with the Standard Model 



5 



energy limit. The Standard Model including a Higgs sector is precisely renormalizable, 
but if the Higgs is decoupled from the theory by taking its mass to infinity, Mh — > oo, 
we should recover the divergences of the Higgsless case. One can therefore think of Mh 
as a soft effective cutoff corresponding to the scale of new physics [20] . Any large Higgs 
contributions must appear experimentally as deviations from the tree-level expectation, 
which can be summarized in terms of Veltman's "p-parameter." In principle, this could 
receive corrections from the Higgs sector of the form 

p= = l + aqq^— l + aiQ^ n— 1 + (1) 

^ M| cos2 9 °^ M| M| ' ^ ^ 

where oq and ai are pure numbers which must be calculated, g is a coupling constant, and 
'• ■ ■' denotes the effect of higher-order radiative corrections which we have neglected. The 
current experimental constraint is p = 1.0004;1^q;qqq4 [21], so if 7^ one would obtain 
extremely stringent constraints on Mh- Unfortunately, in the Standard Model it turns 
out that oq = [221 123], leading to a considerably weaker bound Mh ^ 215 GeV. This 
effect occurs in all Standard Model observables and has become known as the screening 
theorem, because it protects low-energy observations from the effect of coupling to a large 
phase space of scalar Higgs states. It has been shown that the screening phenomenon 
extends to all orders in the loop expansion in the limit Mh 00 [211 [25], [19] . 

The same principles apply to any light scalar field. What happens if Standard Model 
particles are permitted to radiate into dark energy states? In the laboratory environment 
where the and Z masses can be measured, the dark energy quanta are typically 
light. In this case, we must expect contributions to electroweak observables of the form 
described by Eq. ([1]), with the Higgs mass Mh replaced by whatever scale M determines 
the size of the phase space of available states, and the coupling g'^ replaced by whatever 
quantity sets the interaction strength of dark energy with ordinary matter, which is 
typically a number of order M|/M^. It then becomes extremely significant whether 
dark energy exhibits a similar screening effect, for if oq 7^ then p will generically 
receive corrections of 0(1). Such large corrections could easily lead to an unacceptable 
conflict with precision electroweak data. On the other hand, if the dark energy does 
exhibit screening then the corrections to p are roughly of order 0[{Mz/My lnM^/M|] 
and are therefore very small for any phenomenologically reasonable choice of M. 

We would like to emphasize that there is no reason of principle for the Higgs or 
any other scalar species to exhibit this sort of radiative screening. In the Higgs sector, 
a so-called "custodial" global SU{2) symmetry becomes exact in the limit where the 
hypercharge gauge coupling gi vanishes [26], which guarantees equality of the vector 
boson masses, but does not guarantee screening [2D]|jj In the absence of any specific 
reason to think otherwise, one must imagine that a generic scalar field theory interacting 

f In their proof that the Higgs exhibits radiative screening to aU orders in the loop expansion, Einhorn & 
Wudka made essential use of the SU (2) custodial symmetry [25]. However, although the existence of this 
symmetry is necessary, it is not sufhcient. An integral part of of Einhorn & Wudka's argument consists 
of a power-counting procedure entirely unconnected with the custodial symmetry, which determines 
where the leading divergences can appear as Mh 00. 
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with the electroweak sector would contribute to Eq. ([T]) with 7^ 0. Although it may 
be possible to fine-tune a model of this type to be consistent with precision electroweak 
observations, this solution would be highly unattractive. Indeed, one would have traded 
an unappealing fine-tuning in the cosmological constant for a fine-tuning in the scalar 
model intended to replace it, and little would have been gained. 

2.2. The interaction Lagrangian 

We will choose to work with a theory of the broken phase of the electroweak force in 
which the photon and the massive vector bosons and Z interact with a single dark 
energy scalar x according to the action 

S=-\j d'a; |25(/3x)(5"W^+^ - d''W+''){daW^ - d,W-) + Aml.BH{(3HX)W+''W- 
+ B{PxWZ' - d'Z^){daZ, - di,Za) + 2mlBH{PHX)Z''Za 
+ B{(3xWA'' - d''A'^){daA, - (2) 

where and Za are the gauge fields associated with the and Z, respectively, and 
Aa is the gauge field associated with the photon. Eq. ^ should be thought of as an 
effective Lagrangian valid after integrating out the Goldstone modes of the Higgs, as 
emphasized by Burgess & London [271 [28] following earlier work in Refs. [291130]. Only 
invariance under the electromagnetic f/(l) gauge group is required. 

The quantities mw and mz are the Lagrangian parameters corresponding to the 
mass of the and Z, which are related via a renormalization prescription to the 
physical masses M\y and Mz- In addition, we have introduced two arbitrary functions 
B{Px) and Bh{Phx) which describe how the scalar x couples to the gauge boson 
kinetic and mass terms. These couplings are associated with mass scales M = l3~^ 
and Mh = (not necessarily identicaljl) which control the relative strength of the 
interaction between dark energy and the weak gauge bosons, and between dark energy 
and the Higgs field respectively. 

Throughout this paper, we assume that the dark energy quanta x have some fixed 
mass M^, which is not subject to renormalization. This is tantamount to treating the 
entire scalar sector as an effective field theory, in which quantum effects have already 
been included, and for which we only wish to assess the influence of radiative corrections 
on the bare electroweak sector. This is appropriate for a phenomenological model such 
as a chameleon, which need not be a fundamental particle in its own right, but rather 
may represent the collective effect of degrees of freedom at high energy which have been 
integrated out of the theory. In any such effective fleld theory it is difficult to maintain 
light scalar masses because quantum corrections will typically renormalize these to the 
scale of the cutoff unless they are protected by a symmetry. This difficulty afflicts all 

I Note that Mh is not the Higgs mass, which was discussed in §2.1l but does not appear in the remainder 
of this paper. 
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Figure 1. Contributions to the decay width of the neutral Z boson. In (a), an on- 
shell Z decays to a fermion-antifermion pair //. In (b), the decay is precipitated by 
emission of a dark energy particle, x, leaving the original Z in an off-shell excited state 
which subsequently decays to //. If the final-state x is not observed, these processes 
cannot be distinguished and therefore both contribute to the decay width into //. 



particulate theories of dark energy equally, and we have nothing new to contribute to 
this debate. 

The coupling functions B and Bh are unknown, although they will be subject to 
certain restrictions if we wish the dark energy field to exhibit an acceptable chameleon 
phenomenology. We will not impose any such restrictions, except to observe that the 
coupling functions for the W^, Z and 7 kinetic terms must be the same if Eq. ([2]) is to 
descend from an unbroken gauge-invariant theory of SU{2) x f/(l) at higher energies. 
Moreover, the coupling functions multiplying the mass terms must be the same if we 
suppose that the and Z obtain their masses via spontaneous symmetry breaking, 
and that the Higgs sector consists of a minimal SU{2) doublet. Since we wish to retain 
both these phenomenological successes of the Standard Model, we are left with at most 
two free coupling functions. In many cases, however, we expect that Eq. ([2]) will not 
have a UV completion unless these couplings are the same, because the longitudinal 
polarizations of the Z and are associated with Goldstone modes of the Higgs. 

3. Electroweak precision observables 

3.1. Constraints from Z decay 

Let us first consider corrections where some dark energy quanta are present in the final 
state. These corrections can be considered as a form of "dark energy bremsstrahlung" . 
Since the final-state dark energy particles escape the detector and are not observed, such 
reactions look like extra contributions to the cross-section for the corresponding bare 
Standard Model process. Among the best-measured of these is the width for Z decay, 
depicted for decay into a fermion-antifermion pair // of common mass Mf with and 
without dark energy dressing in Figs. [T](a) and (b) respectively. In the dressed process 
(b) the on-shell 4-momentum of one outgoing fermion (which we label '2' by convention) 
is fixed by conservation of 3-momentum. The energy of the other fermion is determined 
by energy conservation in terms of a quadratic equation to be given below. We show in 
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Appendix A that the differential contribution to the Z decay width from emission of a 
single dark energy particle of energy into a solid angle dfi-^ satisfies 



M 



xff 



where M = j3~^ is the dark energy coupling scale, and x and y are defined by 



(3) 



X = 



y 



Ml 



The outgoing dark energy scalar is taken to have 3-momentum q. 
dimensionless "hatted" energies and momentum according to the rules 
Ei 



Ei 



q 



Mz 

q 



(4) 
(5) 

We introduce 

(6) 
(7) 



where i G {x, 1,2}. The quantity measures the degree to which the intermediate Z* 
is off-shell, and satisfies 

r' = -l + 24 - y\ (8) 

It is equal to —1 for an intermediate Z which is precisely on-shell, although in this limit 
the finite width of the Z cannot be ignored. The energy Ei must be a solution of the 
quadratic equation 



Ef [cos'' 9{El 



Ei(i-4)(i + i/2_24 



y ) 

l{l + y^-2E^f + x\cosH){El-y') 



(9) 



where 9 is the angle between q and the 3-momentum of fermion 1. Although two 
solutions for Ei exist, one is always spurious. The solutions change roles at 9 = tt/2. 
Moreover, J is a Jacobian arising from fixing Ei to be a solution of Eq. ([9]). It is defined 
by 



J 



1 + ^1 



l + (4 



y 



2\l/2 



(4 _a;2)-l/2 COS 9 



- El 



The matrix element M.ff satisfies 



Mfj = VT^^{QgL9Rx' + {gl + gl){l - x')} 



(10) 



:ii) 



where gL and gn are the left- and right-handed couplings of the fermion species to the 
Z] and ■M.y.ff is a complicated function whose form is determined in Appendix A and 
which can be read off from Eqs. ( 1A.9[) -( IA.10I) or Eqs. ( ]A.14l) -( lA.15[) . After integrating 
over 4 3.nd the solid angle Eq. ((31) determines the cross-section for any dressed 
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process in terms of the bare standard model cross-section. In a generic model without 
fine-tuning, for which i?' ~ 1, this rate takes the form 

T{Z-.xff) _ 1 M| 

and I^fj is found to be numerically of order I^jj ^ 0.2 for a wide range of fermion 
masses and couplings. The width of the Z into visible particles is predicted to be 
Tz = 2.4952 GeV within the Standard Model, with a small theoretical error. Its 
measured value is Tz = (2.4952 ± 0.0023) GeV [21], implying that any enhancement 
due to dark energy will be compatible with observation only if M > 0.66M^ ~ 60 GeV. 
Moreover, our neglect of the Z width means that this is a conservative over-estimate. 
Thus, under the very mild constraint M > Mz it seems clear that there will be no 
disagreement with the data. Processes similar to Fig. [11(b), but with emission of more 
than one dark energy particle into the final state, are suppressed by extra powers of 
{Mz/M)\27c)~\ 

Dark energy bremsstrahlung could have consequences beyond enhancements to 
decay widths and cross- sect ions of the sort calculated above. Soft bremsstrahlung effects 
could be significant in QCD if they initiated jet formation by destabilizing quarks 
or gluons, or if their aggregate effect could be resolved by partons participating in a 
sufficiently hard collision. However, such effects are likely to be important only if the 
dark energy couples at a very low scale. We can estimate that the S-matrix element 
for any bremsstrahlung event should controlled by the square of the single-chameleon 
coupling constant, of order Mj/M for a fermion of mass Mj, and a phase space factor 
of order In s/M'^, where s ~ M^-y^ is the usual Mandelstam variable and ^ 10~* eV 
is the dark energy mass in the beam pipe [311 ES] . The logarithm is roughly of order 
10^. A significant effect can occur if the product {Mf/M)'^\nM^y^/M^ ~ 1, but unless 
the dark energy scalar is very light this combination is generally negligible whenever the 
coupling scale M is modestly larger than the mass of the fermion species in question, 
of order M > lO^M/. 

3.2. Oblique corrections 

In addition to bremsstrahlung processes, the perturbation theory constructed from 
Eq. ([2]) describes processes by which Standard Model particles may radiate into an 
intermediate state containing an arbitrary number of dark energy quanta. If we exclude 
reactions in which dark energy particles are present in the initial or final state then 
all such processes are built out of interactions which are already present in the bare 
Standard Model. To study them we should begin with a given Standard Model reaction, 
exemplified for the case of 2 — 2' scattering of light fermions in Fig. [2(a), and account 
for the effect of dark energy activity. This activity can naturally be divided into three 
categories, corresponding to Figs. [2](b)-(d). 

In Fig. [2](b), dark energy loops dress each vertex in the bare reaction with so- 
called daisies, whereas in Fig. [2](c) dark energy quanta bridge between two different 
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(a) 




Figure 2. Classes of dark energy diagrams associated witlr Standard Model reactions, 
exemplified in the case of 2 ^ 2' fermion scattering. Solid lines with arrows represent 
fermions; wavy lines represent the gauge bosons of the electroweak force; and plain lines 
represent dark energy particles. The bare Standard Model process is given in (a). In 
(b), the vertices of the reaction are dressed by daisies which begin and end at the same 
vertex. In (c), dark energy quanta bridge between two diiferent vertices. Corrections 
such as (b)-(c) which depend on the process under study (in this case, depending 
on the initial and final fermion species, and the identity of the exchanged boson) 
are called straight. On the other hand, corrections such as (d) which are universal 
for all processes involving the exchange of a given species of vector boson are called 
oblique. (In principle there are also oblique corrections to the fermion species, but 
typically these do not contribute significantly to observable quantities.) In general, 
the dark energy correction to (a) consists of summing over all possible combinations 
of processes similar to (b)-(d). 



vertices. More complicated bridges, including internal vertices which may themselves 
be dressed by daisies, can also b; constructedfl Together with processes where one or 
more dark energy quantum appears in the final state, these are examples of so-called 
straight corrections which depend on the process which under consideration [33] • 

In contradistinction. Fig. [2{d) represents an example of an oblique correction, 
which involves intermediate dark energy states only in the interior of a gauge boson 
propagator. Once an oblique correction has been calculated for a given species of gauge 
field, it is universal for all processes involving exchange of that boson. In principle, 
these corrections are all equally important and for a general momentum transfer q 



it is a complicated process to compute them. However, we argue in Appendix B 



that the daisies and bridges which constitute the straight corrections are momentum- 



t Note, however, that we do not include loops in which the chameleon interacts with itself: as has been 
said, these are assumed to have been absorbed in the parameters of the dark energy model. 
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Figure 3. Processes contributing to the self-energy of the intermediate vector bosons 
7, and Z. An initial vector boson state, represented by a wavy line, radiates into 
scalar quanta x (represented by a solid line) which are eventually re-absorbed to yield 
a final state characterized by the same quantum numbers and momentum as the initial 
state. 



independent up to terms of order g^/M^, where M > Mew is a dark energy coupling 
scale characteristic of the fermion species which participate. Provided they are the same 
for all species, such momentum-independent terms can be absorbed in a renormalization 
of the Fermi constant, Gp, and are therefore unobservable. We have seen in §3.11 that 
in any phenomeno logically acceptable scenario we expect M ^ M^w, implying that 
the remaining contributions can be neglected in comparison with that of the oblique 
correction. Fig. [2][d), which is present at order g^/Ml^v- Oblique corrections will 
therefore give the most stringent constraints if they turn out to require M > 100 GeV. 

The effect of physics beyond the Standard Model has been studied by many 
authors, and is frequently dominated by oblique corrections. Peskin & Takeuchi [33l [31] 
introduced a simple parametrization of them in terms of three quantities S, T and 
U which quantify the magnitude of corrections near zero momentum transferfl but 
assumed that whatever new physics was responsible for modifying the properties of the 
gauge bosons was heavy. This assumption was later removed by Maksymyk, London 
& Burgess [371 EH], who introduced new parameters V, W and X to quantify the 
significance of radiative corrections around the Z resonance|§| In the remainder of this 
section, we briefly review the parametrization of oblique corrections in terms of S, T, 
U, V, W and X. 

The one-loop obliquely-corrected vector boson propagators are obtained by 
summing over an arbitrary number of insertions of the one-loop diagrams of Fig. [3] 
in the tree-level propagator. In unitarity gauge, where the three would-be Goldstone 
modes supplied by the Higgs doublet have been absorbed as longitudinal polarizations 
of the and Z, the tree-level propagator for each massive vector boson can be written 



where we have defined k = ki = —k2 and the quantum field Xa is built out of the 
creation and annihilation operators corresponding to a vector boson of species A and 
mass m^. The tree-level propagator function satisfies A~^(fc^) = k"^ + m\. The photon 
propagator can be written in an analogous form, with i— > in the function A(/c^) 

I An alternative parametrization was proposed simultaneously by Altarelli & Barbieri |351 136] . 
§ See also Refs. jMjliQj. 



{X%{k,)X^i{k,)) 



i{2'n)H{k^ + k2)5AB Tl"' + 
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and m\ i-^ —k"^ in the tensor prefactor. 

We define the sum of the one-particle-irreducible diagrams which connect an initial- 
state vector boson of species A with a final-state vector boson of species B and carrying 
momentum k to be in^''^(A;^)/(27r)^. Since the Z and 7 are electrically neutral they 
are permitted to mix beyond tree-level, which would correspond to a non-zero vacuum 
polarization H^''^. However, inspection of the interactions in Eq. ([2]) shows that Eq. ([2]) 
does not induce extra mixing and we can set 11^''^ = 0. 

With this simplification, the full propagator can be resummed using the Schwinger- 
Dyson equations. The result is that the propagator function A in Eq. (fT3!) should be 
replaced by a resummed function A', which for each species A satisfies 

A'(A;') = , (14) 

where we have written 

= v^'nfsik') + k^k'n^iUk'), (15) 

and, for external states which consist only of light fermions of invariant mass-squared 
M|, Eq. (fT4|) is valid up to corrections of order Mj/M^ which we neglect. Therefore, 
the quadratic term will not appear in the remainder of this paper, and to simplify 
notation we write all subsequent formulae in terms of the abbreviation Hab = ^^ab- 

3.3. The S , T , U, V and W parameters 

In the absence of radiative corrections, the Standard Model entails the existence of 
simple relationships among the observables of the theory. Since there are three free 
quantities which parametrize the broken phase — the two gauge couplings gi and g2, 
together with the Higgs vacuum expectation value — it is necessary to take three masses 
or couplings from experiment. Once this so-called 'input parameter set' has been 
selected, all other observables can be expressed in terms of the chosen three. In the 
electroweak sector it is conventional to choose the input parameter set to comprise the 
fine structure constant a, the Fermi coupling Gp, and the Z mass, Mz, which are 
presently the best measured electroweak quantities. 

With the inclusion of radiative corrections, the original simple relationships among 
observables are modified. Indeed, in order to match the precision with which accelerator 
experiments can measure electroweak parameters, it is usually necessary to include 
several orders of radiative corrections which arise purely within the Standard Model. 
It may happen that these corrections are insufficient to account for the deviation of all 
observables from their tree-level values. The remainder must be ascribed to new physics: 
it is only this contribution from new physics which we wish to attribute to the effect of 
a dark energy scalar species. The relevant observables other than {a,GF, Mz} are the 
mixing angle, 6'w, and the mass, Mw, together with any cross-sections or decay 
rates which can be written in terms of all these quantities. At tree level, 6'w and Mw 
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are related to the input parameter set via the rules 

s,„=.w(l-si„=«w)^^j^^. (16) 

Mj. = M|cos2 9w, (IT) 

where cos6'w in Eq- ( jl^ is to be computed from the solution to Eq. ( fT6l) . 

The physical mass of a single-particle state corresponding to a vector boson is given 
by the pole of Eq. (fT4l) . which renormalizes the Lagrangian parameter m^. Therefore, 
the physical mass M4 satisfies 

Ml = Mi (1 - . (18) 

where we have introduced a useful notation in which a tilde, as in Ma, denotes the value 
taken by a quantity in the Standard Model without oblique corrections. At tree-level, 
M\ is simply equal to m\, but Eq. (fTSll continues to apply to leading order in radiative 
corrections even if we allow the vector boson masses to receive renormalizations from 
loops purely within the Standard Model. On the other hand, the Fermi constant Gi;' is 
defined as the coupling of the charged-current interaction at zero momentum transfer 
and receives an obhque correction [331 El] 

Likewise, the fine structure constant measures the electromagnetic interaction at zero 
momentum transfer and receives an oblique correction from the photon self-energy. 



a = ai^l + n^T,(0) j , (20) 

where 11^^ (/c^) = Ilj^{k'^)/k'^. Eqs. (|T9|) and fl20l) apply even if we allow Gp and a to 
receive corrections from pure Standard Model loops. It follows that we can write 

and 

where we have introduced the useful abbreviations sw = sin6'w and cw = cos6'w, and 
the parameters S, T and U are defined by [33l [3ll l37j ||l 

II Certain terms in these expressions change sign depending on the choice of signature for the metric. 
Under reverse of sign convention (which gives the timehke convention widely used in particle physics, 
in comparison with the spacelike convention adopted in this paper), the formulae for S, T and U 
should be modified by reversing the sign of each mass-square M^, together with extra signs for each 
explicit factor of or dfc^. This explains the difference in signs between Eqs. ([23|) -l[25 |) and the original 
references, which used the signature (+,—,—,—). Note also that in theories where the Z and 7 mix 
beyond tree- level, S and U receive extra contributions. For details, see Refs. [33 l [34 l [37]. 
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^ nzz(o) _ nw(o) .... 
« , m - nw(o) - niyiy(-M2,) 



^ , ([/ + 5) = ::r'' - n,,(o). (25) 

Experimentally observable quantities such as the Veltman p-parameter, Eq. ([T]), can be 
written in terms of S, T and U. 

Electroweak data is not limited to measurements of the and Z masses and the 
mixing angle, but includes cross-sections and decay rates. The standard LSZ formula 
|5T] implies that the first-order shifts from oblique corrections in these quantities can be 
obtained from their tree-level values together with appropriate multiplication by wave 
function renormalization factors Z^, defined for each species of massive boson A by the 
rule 



(26) 



To take account of these factors, it is necessary to introduce two further parameters V 
and EH Em SoT 



2^ 



nzz(O) - Uzz{-_Mll^ ^27) 



Oblique dark energy corrections to all purely electroweak observables can be written in 
terms of S, T, U, V and W. 

These parameters have simple physical interpretations. S" is a measure of the 
difference between the wavefunction renormalization of the Z boson and the photon, 7. 
In an interacting theory, a state prepared with definite particle content and momentum 
at some early time may not manifest the same content when probed at a later time 
because the particles may radiate into any other states to which they couple. The 
probability for this to occur is quantified by the wavefunction renormalization. 

T is a measure of the extra isospin breaking at zero momentum which is contributed 
by new physics. This difference manifests itself in the relative strength of the charged- 
and neutral-current interactions. The precise balance between these interactions may be 
upset by coupling to the dark energy scalar, but in the Standard Model with a minimal 
Higgs sector T is unlikely to receive large corrections unless isospin symmetry is broken 
explicitly at tree level. Similarly, f/ is a measure of the difference between the and 
Z wavefunction renormalizations. Finally, V and W quantify the difference between the 
wavefunction renormalizations of the Z and bosons, respectively, on the mass-shell, 
compared with zero momentum. In what follows, we will see this structure emerge 
explicitly from our analysis. 



% In theories where the Z and 7 can mix beyond tree-level, it is necessary to introduce a third new 
parameter, X. See Ref. [23 . 
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4. Vector boson vacuum polarizations 

To evaluate the S, T, U, V and W parameters, one requires an explicit expression for the 
vector boson vacuum polarizations. In this section, we obtain the necessary self-energies 
by calculating the two one-loop diagrams in Fig. [3l 



4-1- Feynman rules 

Our detailed information concerning the properties of the and Z bosons comes 
mostly from the LEPII experiment, which created these particles abundantly in head- 
on e"'"e~ collisions. The W^s and Zs synthesized in this way were produced at rest in the 
beam-pipe and spent their entire lifetime within its vacuum, before decaying into other 
particles which could subsequently be detected. In the environment of the beam-pipe, 
we can assume that the scalar field has a constant vacuum expectation value x together 
with small excitations 5x- To obtain the one-loop vacuum polarization, it is necessary 
to describe the interactions of the and Z to order 5%^. For interactions involving a 
and W~ the relevant vertices are: 



k 



3 



^B'f3[r,'^\k2-k^-Ww)-klkl\, (29) 



"a 




B"f3 



- [v''\h ■ ks - eml,) - k\k%\ , (30) 



where B' = B\l3x), B" = B"{Px) together with equivalent definitions for Bh; the 
spacetime inner product is denoted p ■ q = p'^Qa for any two four- vectors p°' and q"", and 
we have defined quantities 7 and e according to the rules 

7=^% (31) 
' B' p ^ ^ 



With this choice of Feynman rules, the diagram of Fig.[3](a) corresponds to a vacuum 
polarization of the form 



/?2 B'^ ^ dK K 



^ww{.k ) = / dx 



2 



-{2k' + YM^) + {xk' + ^M^) 



(33) 



87r2 B Jo Jo (k2 + E2)2 

where a; is a Feynman parameter, and we have Wick rotated to Euclidean signature 
before replacing the Euclidean volume element by 27r'^K^ dn. The momentum scale A 
is a sharp cutoff which regulates the maximum Euclidean momentum permitted to 
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circulate in the loop, and therefore determines the size of the phase space of scalar 
states to which each W'* couplesH Finally, S is an abbreviation for the quantity 



= x{l - x)P + (1 



(34) 



In writing Eqs. (!33|) - (!34|) we have freely replaced by M^, since the correction this 
induces is formally of higher order in the loop expansion. 

The diagram of Fig. [3](b) gives a somewhat simpler contribution. 



(3^ B 



12 



87r2 B 



dx 



K' d/t 
k2 + M2 2" 

A 



'w)-> 



(35) 



where f2 is a dimensionless combination which measures the curvature of the coupling 
function B in the vacuum, 

B"B 



B 



12 



(36) 



We also require the vacuum polarization for the Z boson and the photon, 7. 
However, no further calculation is required since the relevant Feynman rules can 
be obtained from (I29 l) -( l30l) . and the necessary vacuum polarizations can likewise be 
obtained from Eqs. (!33|) - (!35!) . Since the 7 and Z are their own antiparticles, each 
vertex in (129|) - (!30|) acquires a symmetry factor of 1/2. To obtain the correct vacuum 
polarizations, one makes the replacement Mw ^ Mz in Eqs. ( I33i) - (l35l) for the Z, and 
Mw for the photon. 

Assembling these terms and carrying out the k integrals, it follows that the vacuum 
polarization for each species of boson satisfies 



/2 



dx 



B 

+ {xe + jMl)^ 



2fc2 + 72MI 



A2 + 



A2 A2 

YA2 + S2 



S^n 1 



A2 

S2 



2A2 



A2 1 , 

Ts^ + 2^" 



A2 



A2 
2 



— In 



1 + 



A^ 

M2 



(37) 



4-2. Effective Lagrangians for the vacuum polarization 

Eq. (1371) is a complicated expression from which it is difficult to extract the important 
qualitative features of the oblique corrections. To do better, one can analyze 11^^ (fc^) in 
terms of an effective Lagrangian which would give rise to the same vacuum polarization. 

f Power-law divergences in A, if they exist, are likely to violate gauge invariant although logarithmic 
divergences should be physically meaningful. Also, loop calculations in unitarity gauge are known 
to overestimate power law divergences in certain circumstances. These issues were addressed in 
Refs. [27l [28]. In the present case it will turn out that we require only the logarithmic terms. If 
any power-law divergences were present, however, then it would not be possible to interpret the result 
as a quantitative prediction. Instead — provided such powers were compatible with naive dimensional 
analysis (which excludes the possibility of overestimation) and the gauge symmetries of the model — the 
correct interpretation would be that the calculation under discussion was sensitive to the details of UV 
physics. 
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A. Low energy, massive vector bosons. Consider first the limit <^ M, 

each species of massive vector boson A one can make the expansions 



2 For 



A' 



A2 + S2 



oo 



n=l 



x(l — x) 



and 
In 



A2 



In 



cr 



a" 



^_-|^'jm.(n+l)+l 



m=l 71=1 



mcT 



2m 



A2 



a;(l — x)k'^ 



x)Ml + xMl 



where we have defined cr^ by the rule 

l + (l-a; 



a 



Ml Ml 
A2 ^''a^' 



(3^ 



(39) 



(40) 



In particular, o"^ ~ 1 whenever the scale of the cutoff, A, is much larger than the 
electroweak scale Mew ~ ^a- Eq- ( |38l) is an expansion in powers of /c^/A^. In an 
effective field theory, these contributions would come from a tower of non-renormalizable 
operators suppressed by the cutoff scale, although one should remember that whenever 



these operators become important the bridge corrections discussed in Appendix B 



will also make a significant contribution. On the other hand, Eq. (139|1 amounts 
to an expansion in powers of /c^/M^ 1^1 These contributions would come from non- 
renormalizable operators suppressed only by the electroweak scale. As we increase 
the momentum which is transferred through the gauge boson propagator from zero, 
we expect to see corrections enter at the scale |/c^| ~ M^, followed by another set of 
corrections at the cutoff. 

Collecting these expressions one finds an expansion for 'n.AA{k'^)i which yields 



3_ 

M2 



Mlao + aaA;^ + a^k'^ + O 



k 



,2 



« M|w) 



where g is an effective dimensionless coupling constant defined by 

1 fi'2 



9 



87r2 B 



(41) 



(42) 



the mass scale M is M = (3 
72 \ 72^ 



ao = 


A2 / 

^ V 


a2 = 


4 ^ 


0^4 = 


12 



+ 



Ml 



16 



as before, and the coefficients a^, for i G {0, 2, 4}, satisfy 
" A2 



6 In 



Ml 



A^ 

Ml 



144 



-1 + 



A2 



Ml 



6[7(12 - 7) - 6] In ^ + 7(36 - 57) - 18 + O ( 



Ml 



Ml 



(43) 
(44) 
(45) 



X The series expansion in Eq. (|39[) can be integrated term-by-term in x, producing an expansion in 
powers of k"^ /Ml with coefficients which involve hypergeometric functions of M^/M^. When expanded 
in powers of this ratio it is possible that logarithms of Ml/ Ml are generated, although suppressed by a 
positive power of M^/M^. It follows that the scale at which this tower of non-renormalizable operators 
becomes significant genuinely is around the electroweak scale, |fc^| ~ Ml. 
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We could equally well have obtained this vacuum polarization if we had started from 
an action of the form 



+ corrections at Mp 



EW 



(46) 



and calculated only to tree level, where ip represents any polarization of the vector boson 
of species A, and the corrections at M^-y^ take the form of a tower of non-renormalizable 
terms suppressed by powers of Mew Note the unsuppressed non-renormalizable term of 
the form ipd'^yD, which is symptomatic of the fact that our starting Lagrangian, Eq. ([2]), 
did not describe a renormalizable quantum field theory. 

A good deal of information can be obtained from inspection of the effective 
action P6|) . The relevant operators are the kinetic term cpd'^ip and the mass term 99^, 
which both receive corrections quadratic in the cutoff A. The mass is prevented from 
receiving corrections which scale faster than A because gauge invariance is restored when 
Ma —>■ 0, and in this limit the mass must not receive quantum corrections so that the 
Ward identity is preserved. Indeed, it follows from Eq. fHT]) that the T parameter can 
be written 

aT" = ^ (ao,z - ao,w) ' (47) 

and therefore that all quadratic divergences cancel in this quantity. It is clear from 
Eq. (jlHD that this cancellation is a direct consequence of the restoration of gauge 
invariance in the limit Mz, My/ 0. 

B. Low energy, massless vector bosons. A similar procedure can be applied to 
find an effective Lagrangian for the photon self-energy in the low-energy limit. The 
vacuum polarization is obtained from Eq. fl?r|) after the replacement Mw \—>- 0, after 
which the expansions (l38!) - (!39l) continue to apply, with cr^ substituted by the alternative 
combination r^, which satisfies 

r'^l+x^. (48) 

However, the roles of these non-renormalizable operators are subtly changed. Eq. fl38l) 
can still be interpreted as a tower of corrections at the cutoff (which we again caution 
will be accompanied by significant bridge contributions), but Eq. (139|) now represents 
corrections at the scale of the dark energy mass, \k'^\ ~ M^. If we discard these 
corrections, it follows that the effective Lagrangian we obtain will be valid only in the 
limit <^ M^. Fortunately, for finite this is sufficient for the purpose of obtaining 
the oblique parameter S. 
In this limit, one finds 



n^^(A;2) - 



2 

(l^'l <^v) (49) 



X 
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where the coefficients 62 and ^4 satisfy 



A2 

62 = — (l-fi) + 0(M, 



5, = — 



6 



O 



V A2 



(50) 
(51) 



Within its range of vahdity, this expansion can be interpreted in terms of the effective 
Lagrangian (l46l) . In particular, note that (as expected), no mass term is generated 
owing to gauge invariance. 

C. Energies near the resonance, massive vector bosons. To obtain we 
require information about H-AAik'^) in the region where it approaches the resonance 
at k"^ = —M\. This can be studied by setting fc^ = — M^ + g^, and making an expansion 
in powers of /M\. When expanded in this way, it is less straightforward to interpret 
n(A;^) as an effective Lagrangian. However, some of our understanding concerning the 
meaning of each term can be carried over. 

Eqs. f l38l) - (!39l) . giving expansions in terms of non-renormahzable operators, 
continue to apply with the replacement t— > a^, where for each species A of massive 



vector boson we have defined 
= 1 + (1 



We find 



M\ao + a2q^ + a^q^ + 0(^ 



Ml) 



(52) 



(53) 



where the coefficients dj, for i G {0,2,4}, satisfy 



A^ 



A2 

a2 = — 



7 



-n{e-l)-l 



Ml 



/ A^ / 

+ (^3[7(57 - 6) + 4] In ^ + 4(7(47 - 9) + 5] + O ( 

;i - fi) + ^ (^3[7(ll7 - 16) + 4] In ^ + 7(677 - 128) + 59 + O ( 



(54) 



M 



1 , A2 
"^=12^"m! 



Ml 



(55) 
(56) 



It is now possible to give expressions for the remaining oblique parameters S, V 
and W in terms of these effective quantities 



aS 



4^2 p2 



(ao,z - «o,z - ^2) 



aV -- 
aW 



g_ 

M2 



— (d2,z + ao,z - tto.z) 



3 , - \ 



(57) 

(58) 
(59) 
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where we have dropped contributions from the non-renormahzable operator <fd^(f since 
these never lead to quadratic divergences. It is now clear from inspection of Eqs. fl57|) - 
( 1591) together with Eqs. (03]) -(05]), (l50])-(l5ll) and (I55])-(I56D that all quadratic divergences 
cancel in S, T, U V and W. 

5. Discussion 

5.1. When are quantum corrections screened? 

This cancellation is not an accident, but is partly a consequence of gauge invariance and 
partly depends on the structure of gauge boson-lepton couplings within the Standard 
Model. 

The available phase space which sets the size of the loop correction is determined 
by the couplings {B, Bh} and the mass of the boson, which is an infra-red effect. Most 
of the phase space volume will be concentrated near the ultra-violet region, in spherical 
shells of large Euclidean four-momentum. Coupling to these shells corresponds to a 
process where a propagating intermediate vector boson radiates into a hard chameleon 
and boson pair. From the point of view of this pair, the original vector boson behaves 
as if it were massless, and the effect of mass splittings between W^, Z and 7 becomes 
irrelevant. Therefore, because gauge invariance requires that W^, Z and 7 couple to 
the dark energy in the same way at zero mass, we expect no difference in the manner 
in which any of these gauge bosons radiate into momentum shells at Euclidean four- 
momenta which are large compared with Mz- 

Assuming our choice of input parameters, this is sufficient to screen all 0(1) effects 
in contact interactions of a single electroweak gauge boson with exactly two fermions — 
which is the only type of interaction which occurs in the electroweak sector, excluding 
interactions with the Higgs. The input parameters were chosen to be the Z mass, 
Mz, together with the fine structure constant, a, and the Fermi constant, Gp, which 
measure the strength of the electromagnetic and charged- current interactions at zero 
momentum, respectively. Operationally, both a and Gp measure a combination of some 
dimensionless coupling constants and a propagator at zero momentum: for a this is the 
photon propagator, whereas Gp measures the W propagator. The oblique corrections 
can be of two kinds. Firstly, for processes involving a Z particle, the strength of the 
neutral-current coupling is not measured hj Gp but can be obtained from it by a shift 
measured at zero momentum. This is the purpose of the T parameter. Secondly, a 
wavefunction renormalization of gauge boson lines may be necessary, which depends on 
properties of the propagator near Euclidean momentum of order Mz- The difference 
between the wavefunction renormalization of the Z and W propagators evaluated at zero 
momentum and at momenta near Mz is measured by V and W, respectively. Finally, 
S compares the zero-momentum Z and 7 propagators and therefore plays the same role 
for the photon as T does for the Z, while U measures the difference between the 
and Z propagators at zero momentum. 
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Figure 4. Current collider constraints on the coupling scales M and Mh, 
associated with dark energy interactions with the electroweak gauge and Higgs sectors, 
respectively. The interior light green region is compatible with current precision 
electroweak data at la, and extends indefinitely to large AI and Mj^. Also shown 
is the 2a region in darker green. 

All these shifts depend on a comparison of the phase space available to two different 
gauge bosons, or between the same gauge boson at different momenta. As we have 
seen, gauge invariance guarantees that the phase space available to all gauge bosons 
is the same at large Euclidean four-momentum, so differences can only arise from the 
interior shells of momentum space where the mass splitting between the electroweak 
gauge bosons can no longer be neglected. Differences in this region can not lead to 0(1) 
effects if the mass scales {M, Mh} associated with dark energy are much larger than the 
electroweak scale. It follows that large effects from radiative corrections are screened. 
However, this depends essentially on the fact that y/a and y/Gp include one gauge 
boson line, and all other processes subsequently involve vertices which also include only 
a single ingoing or outgoing gauge boson. 

The calculation of the oblique corrections in the previous sections was done only to 
one loop. In principle loop corrections of any order could contribute 0(1) corrections 
but we expect that the screening of oblique scalar field corrections to the gauge boson 
propagators occurs at all orders. 

For a dark energy species which selects its mass via a chameleon mechanism we 
depict the current collider constraints on the mass scales M, Mh in Fig. |H These can 
loosely be summarized as M, Mh ^ 1 TeV, which is stronger than the constraint which 
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Figure 5. Current constraints on the mass, M^, of a dark energy particle wliose 
interaction with ordinary matter is characterized by a scale M > Mew The interior 
light green region is compatible with present data at la, whereas the 2a region is 
shown in darker green. Compare with Fig. 8 of Ref. [H]. 

follows from the decay width of the Z into visible particles. Clearly, neither constraint 
is competitive with present bounds from optical or axion-search observations [121 US] • In 
Fig. [5] we show the same constraints as a function of the dark energy mass, M^, and its 
interaction scale M, without assuming that My is determined by some chameleon-type 
mechanism. 

5.2. Future prospects 

Any future linear collider will measure electroweak precision observables with remarkable 
accuracy [44j, but if weak couplings imply it cannot produce dark energy particles 
directly then the most important discovery mode will come from sensitivity to radiative 
corrections at high energy. For any electroweak processes sensitive to the diverging phase 
space of dark energy states at large Euclidean four-momentum, the discovery reach of the 
ILC would not be limited by the smallness of the coupling unless new physics operating 
at lower energies could quench the contribution of dark energy loops. An example of 
such new physics could be the appearance of a chameleon superpartner at some energy 
MgusY, if ^susY ^ M, where M is the characteristic mass scale with which dark energy 
couples to the gauge sector. On the other hand, if dark energy radiative corrections are 
screened, then contributions to electroweak precision observables fall with the mass 
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Figure 6. Penguin diagram contribution to the decay i?" K^4). 

scale of the coupling like {M^^ / MflxiM^ / Ml^. The most st rinerent constraint on M 
currently derives from the polarization of light from astrophysical sources, which was 
studied in Refs. [151 US] and leads to the lower limit M > 10^ GeV. It is unlikely that 
such small corrections could ever be observed at the ILC. 

Since detection of electroweakly interacting dark energy at e^e~ colliders will be 
challenging, it is natural to consider what can be achieved at hadron-hadron colliders 
such as the LHC or the Tevatron. Although W^s and Zs are produced by such colliders, 
the problem of backgrounds and the difficulty of kinematical reconstruction of the final 
state at a hadron collider mean that constraints from pure electroweak processes are 
likely to be inferior to those from a future linear collider. However, hadron colliders are 
sensitive to other channels in which new physics can appear. One particularly interesting 
window on new physics may be provided by rare decays of B mesons, which are bound 
states of a bottom quark b with some other quark q in the combinations bq or bq. Such 
mesons can decay via fiavour- changing neutral currents which are forbidden at tree-level, 
but give rise to decays such as — > K^(j) (where B^ is the neutral B meson composed 
of an anti-bottom/down pair bd) when loop diagrams such as the so-called "penguin" of 
Fig. [6] are included. Rare processes of this type give a comparatively clean signal of new 
physics. Unfortunately, it does not appear likely that conformally coupled dark energy of 
the kind studied in this paper could manifest itself in this way. At large Euclidean four- 
momentum, where the internal line in Fig. [6] could be expected to receive sizeable 
dark energy modifications, the quarks to which it couples are effectively massless and 
the loop is fiavour-independent. When summed over all quarks which can circulate in 
the loop, the unitarity of the Cabbibo-Kobayashi-Maskawa matrix implies that this 
dominant fiavour-independent contribution suffers an exact cancellation: this is the so- 
called Glashow-Iliopoulos-Maiani mechanism. We can estimate that this mechanism 
allows any dark energy contribution, coupling at a scale M, to contribute at most at 
relative order m^/M^, where ^ 175 GeV is the top mass. 

Is there any way to avoid the screening of large radiative corrections? This can 
only be done if at least one coupling constant measured in a low-energy interaction can 
appear in a different context in some other process. Remarkably, the Standard Model 
does allow for this possibility. If we assume a minimal Higgs sector, there are three- 
and four-body interactions of the massive electroweak gauge bosons with physical Higgs 
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quanta which are described by the action 

S = - J (^^'^^G^H + Gf^^ {2Ml.W+^W- + MlZ^Za) , (60) 

where H is the physical Higgs field. This must itself be subject to oblique corrections 
which only involve the coupling Bh- There is no reason to expect that the shifts 
necessary to bring Gp and the H and gauge boson lines to finite momenta will be 
independent of the ultra-violet region of momentum space. However, these large effects 
are undetectable until the interaction of the Higgs with at least one of the massive 
gauge bosons becomes accessible to experiment. Even when this is possible, the details 
will depend sensitively on the mechanism of electroweak symmetry breaking chosen by 
Nature. For this reason we defer investigation of such processes, although we note that 
in the case of a minimal SU{2) doublet it is possible to verify that one could perhaps 
expect an 0(1) modification of the Higgs production rate via weak boson fusion at the 
LHC. 

The insensitivity of electroweak collider experiments to weakly coupled dark energy 
is frustrating given the inability of cosmological observations to place bounds on this 
region of parameter space. Although the search for astrophysical constraints has been 
fruitful [To], [15], it is difficult to imagine any astrophysical processes which would be 
sensitive to energy densities of order (lO^^GeV)^ or above which were attained only 
during the very early universe. For example, one might have imagined that small 
perturbations imprinted in the dark energy scalar during an epoch of primordial inflation 
would lead to interesting effects in the late universe. Unfortunately, it appears that 
dark energy scalars of chameleon-type generically roll rapidly to their potential minima 
during the first few e-folds of inflation, where they remain for the duration of the 
accelerating era [5]. In the minimum, the dark energy field is heavy and is not excited 
by the infiationary expansion. For this reason, it does not function as an isocurvature 
field and cannot source evolution of the curvature perturbation, which might have led 
to interesting constraints from the spectral index or non-gaussianity. Moreover, the 
curvature perturbation is screened from possible non-perturbative effects because the 
chameleon vacuum expectation value is fixed [16]. On the other hand, if the chameleon 
vev shifted appreciably after inflation, it could potentially amplify the steep blue 
spectrum of perturbations imprinted on the U (1) hypercharge field. If this amplification 
were too dramatic, it would lead to an unacceptable collapse of hypercharge fiuctuations 
into primordial black holes at the end of infiation, in confiict with observation |47j . 
However, in practice the chameleon vev does not change sufficiently for this to provide 
an interesting constraint. 

6. Conclusions 

In this paper, we have studied the prospects for collider physics to detect a scalar dark 
energy species which couples conformally to Standard Model states which are charged 
under the electroweak gauge group SU{2) x f/(l). This is particularly interesting for 
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proposals incorporating a chameleon-type mechanism, in which the dark energy field 
may evade stringent bounds on the presence of light scalar bosons by dynamically 
adjusting its mass to be large in regions of high average density. Any such theory of 
dark energy must certainly couple to the Standard Model, although it is not mandatory 
that the dark energy scalar couples to electroweak states. However, if such couplings 
are present, then in view of the theoretical and experimental cleanliness of electroweak 
physics in comparison with hadron processes, one might expect that they would provide 
the most promising means of detection. 

In the minimal Standard Model with Higgs sector consisting of only a single SU (2) 
doublet, this expectation is wrong. Although coupling to a dark energy scalar in 
principle allows for fractional shifts of 0(1) in precision electroweak observables, we 
have shown that in practice such large corrections are "screened," in direct analogy 
with the screening theorem which prevents similar corrections from heavy Higgs states. 
Screening occurs because a majority of the dark energy corrections are absorbed in 
the input parameters {a, G^, M^}, with only small splittings between the remaining 
observables which arise from the infra-red region of momentum space. The structure 
of the Standard Model also plays an important role, since all relevant vertices involve 
precisely two fermions and a single gauge boson. On the other hand, if it is possible to 
observe Higgs processes at the LHC then we would expect 0(1) corrections to the Higgs 
production cross-section via weak-boson fusion. 

Although we have carried our explicit calculations only to one loop, we expect that 
screening of oblique corrections persists to all orders, any of which could contribute 0(1) 
effects as a matter of principle. This is important in establishing the consistency of dark 
energy theories with existing collider experiments, but also implies that the dark energy 
discovery potential of future e'^e~ colliders such as the proposed International Linear 
Collider may be comparatively limited, unless the Higgs can be detected. 

One might also attempt to probe dark energy couplings via hadron processes, for 
which a promising observable might be the so-callcd flavour- changing neutral ciurcnt 
which mediate rare decays of B mesons. Unfortunately, for such reactions the unitarity 
of the Cabbibo-Kobayashi-Maskawa matrix plays a role similar to that of gauge 
invariance in quenching the contribution from shells of phase space at large Euclidean 
four-momentum, where a significant effect could be expected. Other non-electroweak 
effects such as bremsstrahlung also offer an apparently limited discovery potential. 
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Appendix A. Dark energy corrections to the Z width 

In this Appendix we discuss the possibility of enhancements to the observed decay width 
of the neutral Z boson. Such enhancements arise via the process Z xZ* —>■ x4''4^ 
through which a propagating Z emits a dark energy particle x ^"^^ passes off-shell. 
Eventually the off-shell Z decays to a fermion-antifermion pair Z* ipip, but if the 
weakly-interacting x escapes the detector unseen then this reaction is indistinguishable 
from the direct decay Z iptp. 



Appendix A.l. Direct decay 

Let us first recapitulate the textbook calculation of direct decay [48]. This will allow 
us to express the enhancement from dark energy emission as a fraction of the pure 
Standard Model rate. We suppose that a Z particle decays into a fermion species whose 
quanta are created and destroyed by operators associated with the Dirac fields ip and 
according to an interaction Lagrangian of the form 

d''xi,^''Za{gLL + gRR)ij (A.l) 

where the 7*^ are matrices obeying the Dirac algebra {7", 7*} = 2ri'^^, gi and gn are 
arbitrary left- (respectively, right-) handed couplings, and L and R are projections onto 
the left- (respectively, right-) chirality halves of a spinor in Dirac's representation, 

^^1+75 (A.2) 
2 2 ^ ^ 

We use 75 = — i7°7^7^7^, which has unit square 7I = 1 and commutes with all other 
7-matrices. The projection operators L and R obey = L and R^ = R, together with 
the orthogonality relation LR = RL = 0. We will also use the parity transformation 
operator /3 = 17°, obeying = 1, in terms of which /3(7°')^/5 = —7" and /S'y^P = —75. 

Unpolarized decay proceeds according to the diagram of Fig. [D^a). To obtain the 
overall rate, one averages over the three polarizations of a massive spin-1 particle and 
sums over the two spin states of each outgoing fermion. The differential decay rate per 
unit volume of phase space, dv, available to the final state ipip pair corresponds to 

outgoing ingoing 
spins polarizations 

where ^ k schematically denotes the sum of all ingoing momenta minus all outgoing 
momenta. The Feynman amplitude M^^ depends on the polarization of the initial Z, 
labelled s, and its 3-momentum p, together with the spins of the final-state fermions, 
labelled ai^2, and their 3-momenta ti^2- It is defined by 



[M^.i.]"'" = -77;^ KMtOl'lSii + 9RR)v"{t2)] -2^. (A.4) 
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After performing the spin and polarization sums, yielding a trace over Dirac indices, 
this corresponds to a differential decay rate which can be written 

dr = {2.rs{p - - (4^^ 

X 1 12,,,«M| + {gl + gl) (^^^^^^^ " 2t, ■ t,) } (A.5) 

where ti,2 = (-Ei,2,ti,2) and p = {Ez,p) are 4-momenta corresponding to the out- and 
in-going particles, respectively, and an infix dot '■' denotes contraction in the Minkowski 
metric. All external particles are taken to be on-shell, with 4-momentum conservation 
enforced by 6{p — ti — ^2), and the outgoing fermions each have mass M^. 

Conservation of 3-momentum is sufficient to determine one of the outgoing 
momenta. Moreover, performing the calculation in the Z rest frame, symmetry requires 
that the outgoing fermions have equal energies E12 = Mz/2. In consequence, we 
conclude that the total rate of emission into a solid angle dQ can be written 



dr Mz / Ml ( / Ml\] ^ 

do - Dsif - ^4 r''''4 + + - 4) r ^^-'^ 

Appendix A. 2. Decay accompanied by dark energy emission 

Now we return to the more complicated process where the decaying Z is first pushed 
off-shell via emission of a single x particle and subsequently decays into the observed 
fermion pair. This corresponds to the process of Fig. [1(b). The ZZx interaction vertex 
is determined by (129|) . modified as discussed below Eq. (136|) to obtain the coupling to 
the Z boson. 

As above we label the decaying Z with momentum p and energy Ez, and the 
outgoing fermions with momenta ti^2 and energies -£"1,2. The outgoing x particle is 
taken to have momentum q and energy E^. The total decay rate per unit of phase 
space available to the particles in the final state is given by a formula equivalent to 
Eq. (1A.3I) . with the Feynman amplitude M^^ replaced by a more complicated quantity 
M^^^ which satisfies 

In order to avoid confusion with the parity inversion operator (3 = ij^ we have chosen the 
chameleon coupling scale as M, which elsewhere in this paper has been been synonymous 
with the coupling f3 = M~^. The off-shell interior Z carries 4-momentum r = p — q, 
and G is the 'coupling matrix,' 

Gap = [giL + gRR]^p , (A.8) 
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where {a, (3, ■ ■ ■} label spinor indices. Summing over final-state spins and averaging over 
all three initial-state polarizations, we find 



dr = (27r)^5(p -q-ti-t2] 



1 



^2 + M|)2 (27r)32E^ (27r)32Ei (27r)32E2 6Ez' 



M',,^. (A.9) 



where M.' ^ j satisfies 



PaPd 

Ml 



tr {7^G(-l)(i/2 + M^)/3(G)t(7/)t/5(-i/^ + M^,)} -(A.IO) 



We are adopting the usual Feynman convention in which $ = ^°'Za for any 4- vector Z\ 
and P"c is defined by 



P\ ^ {if\-v ■ r - 7M|) + A"} (^6c + ^) . 



(A.ll) 



The trace over Dirac indices can be evaluated by standard methods. It yields 

trace = r/^^ {^LgRMl - 2{gl + g]^t^ ■ ta} + '^{gl + g\W2t{ + t{t\), (A.12) 

plus a term antisymmetric under the exchange c ^ /, which we omit because it 
disappears after insertion in Eq. (]A.10|) . This trace depends only on the final 
vertex and is common between the direct and accompanied decays. Nevertheless, it will 
not cancel in a ratio between the two, because it depends non-trivially on the Lorentz 
index structure by which it couples to the rest of the diagram. This structure receives 
significant modifications when the Z decay is accompanied by dark energy emission. 
To proceed, we must contract Lorentz indices. We find 



+ Tcrj 



7 - 



p ■ r 

Ml 



1 + 



Ml 



(p ■ r + 7M|) 



+ ^{p-r + lMlf + 



Mi 



p ■ r 

Ml 



{r^ + j^M)Z' 



+ {rcPf+PcTf) 



p ■ r 

Ml 



{r' + I'Ml) 



1 + 



Ml 



{p-r + 7M| 



(A.13) 



This can be contracted with Eq. (]A.12I) for the Dirac trace, yielding a final expression 
for A^'^^^. For convenience of expression, let us write -M'^^^ = A + B. We find 



A 



AMlgLgR-2{gl + gl)t^-t2 



+ r^ 



7 



p ■ r 

Ml 



Ml{r' + YMl) + 4(p ■ r + 7MI 
^2 



2\2 



1 + 



Ml 



{p-r + ^Ml) 



+ 2{p ■ r) 



{p-r + ^Mlf 
^ Ml ^ 



p ■ r 



Ml 



(r' + YMl) 



(p-r + jMl) 



(A.14) 
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and 



B 



2\2 



^i9l + 



= 2{p ■ t,){p ■ t2){r' + YM'z) + 2(ti ■t2)ip-r + ^M'^) 



+ 2(r-ti)(r-t2) 



7 - 



p ■ r 

Ml 



1 + 



{p-r + ^Ml) 



{p-r + ^Mlf 



p ■ r 

Ml 



+ 2[(p-ti)(r 
p ■ r 

Ml' 



t2) + {vt2){r-ti)] 



r' + l'Ml) 



1 + 



Ml 



ip-r + ^Ml) 



(A.15) 



Kinematics. As before, 3-momentum conservation is sufficient to determine the 
momentum of one outgoing particle, which we choose to be t2 without loss of generality. 
Energy conservation determines one further scalar coordinate on phase space, which we 
choose to be Ei. The undetermined part of the 3-body phase space is parametrized by 
the outgoing dark energy momentum q and a pair of polar and azimuthal angles (6', (f)) 
which specify the orientation of ti relative to q. 

Let us obtain Ei as a function of the unconstrained parameters. We work in the 
Z rest frame. Energy conservation requires Ey + Ei + E2 = Mz, and 3-momentum 



conservation fixes t2 = — ti 
equation 



q. Therefore we conclude that Ei must solve the implicit 



E'2 



El + El-Ml + 2^El- Ml^El - Ml cos 6. (A.16) 

To obtain an explicit equation, Eq. (1A.16P can be squared and the resulting relation 
simplified. However, in doing so we introduce an extra solution for which Eq. (1A.16I) 
holds only after the exchange cos 6' 1-^ —cos 6'. The solution is spurious and should 
be eliminated. In practice we will find that the two possible solutions exchange roles 
at ^ = 7r/2, but that the matching is smooth. Following this procedure, the possible 
solutions must solve the quadratic equation 

Ef [4cos2 e{El - Ml) - A{Mz - E^f] + AE,{Mz - E^){Ml + - 2MzE^) 

= {Ml + - 2MzE^f + AMI ^o^' 0{El - M^). (A. 17) 

This can be achieved by substituting the correct solution of Eq. flA.lTp in Eqs. flA.Qp - 
(lA.lOp . However, one must also account for a Jacobian factor associated with 
transformation of the 5-function enforcing energy conservation. To obtain this, define 

S = Ez-E^-E^- E2{E^, E,, 6), (A.18) 

in terms of which energy conservation requires S{£). We can now use a change of 
variables to find 

5[E,-E,{E^,e)] 



\dS/dEi 



(A.19) 
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where Ei{E^,9) is a solution of Eq. (lA.17p . The Jacobian d£/dEi can be determined 
using Eq. (1A.16I) . yielding 

1 + {EI- Mlf/\El - M2)-V2cos^ 



08 




dEi 





1 + ^1- 



Eo 



(A.20) 



where E2 is to be determined by Eq. (lA.16p . In sum, the total unpolarized decay rate 
now satisfies 



dr 



B'H4 



-2 



Ef 



Ml 



d^g 



dEi 



-1 



M' 



(A.21) 



dfii 1927r2(27r)3 (r2 + M|)2 E^E2Mz 

where d^i is the element of solid angle associated with ti. 

To proceed it is convenient to introduce dimensionless small quantities x and y, 
given in Eqs. (jl])-®, which determine and in terms of M^, 

= y/xMz and = ,/yMz. (A. 22) 

Also, we can agree to measure energies in units of Mz, introducing quantities E^ and 
£"1^2 which satisfy Eq. ([6]). Likewise, vectors such as r and q can be rescaled according to 
Eq. ([7]), giving dimensionless vectors r and q. In terms of these dimensionless quantities, 
Eq. (1A.6I) giving the rate of direct decay reads 

dT{Z ipip) Mz 



(A.23) 



where Ai^^ satisfies Eq. (jlTl) with the fermion species / taken to be ip. The Jacobian 
\dS/dEi \ satisfies Eq. f fTOl) . and we will denote it J in what follows. Moreover, -M'^^^ has 
mass dimension [M^], so we can introduce an unprimed quantity A^^V'V' which depends 
only on x, y, the hatted vectors and other dimensionless quantities, and is defined by 



(A.24) 



After these replacements, the rate of unpolarized decay accompanied by dark energy 
emission satisfies 



dr(z ^ x^^) 



1927r2(27r) 



3MI dE^dn^ 



El - x^ 



El 



j(i-E^-Ei)(i + f 



(A.25) 



Taking the ratio of Eqs. (1A.25P and (1A.23P we finally obtain our advertised relation, 
Eq. ([3]), which describes the enhancement due to dark energy emission as a fraction of 
the bare Standard Model rate. As in Eq. ( fT2l) . it is convenient to aggregate that part 
of the enhancement in Eq. ( 1A.25P which is independent of the chameleon coupling M 
into a dimensionless integral This will depend on the mass and couplings of the 

fermion species ip, together with the mass of the dark energy scalar x- particular, 



d^;^ dn^ 



El 



x^ 



El 



M 



IrHm,M,,g,,gn) - j _ _ ^^^^^ ^,2)2 M,, 

We give representative values for /^^^ in Table lAll 



(A.26) 
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fermion species 


massless neutrino 


10 ^ eV neutrino 


511 keV electron 


5 GeV electron 


40 GeV electron 




0.22 


0.22 


0.22 


0.21 


0.007 



Table Al. Enhancement factors for Z decay accompanied by dark energy emission, 
to be interpreted in conjunction with Eqs. p2|) and (|A.26p . As the fermion mass 
increases, the phase space available to any decay products diminishes until it is 
forbidden altogether at the kinematic threshold — Mz/2. The enhancement for 
massless or light particles is very nearly independent of their identity. 

Appendix B. Bridges and daisies: dark-energy corrections to all orders 

Even in the effective field theory interpretation, where loops which are purely internal 
to the dark energy sector are ignored, Eq. ([2]) — together with the assumption that all 
matter fields couple conformally — entails a great many possible corrections to Standard 
Model processes. In this Appendix, we argue that to an acceptable approximation a 
majority of these corrections are subdominant; in this approximation, the leading dark 
energy effect comes from the one-loop oblique correction. This assumption played an 
essential role in determining the dark energy corrections in §3) 

A useful example to keep in mind is the case of the graviton, which also couples 
conformally to matter (and indeed all Standard Model states) with a universal coupling 
function, a/ det {rjab + hab), where rjab is the background metric and hab is the spin-2 
graviton field. This non-linear coupling leads to a network of gravitational bridges, 
daisies and oblique corrections, quite analogous to Fig. [2], which also dress Standard 
Model processes. In the case of gravitons these have very little impact on reactions 
taking place at collider energies; in comparison, the structure of the dark energy 
interactions in Eq. ([2]) allow a small number of diagrams to make an 0(1) contribution. 
Nevertheless, many similarities exist between graviton and dark energy phenomenology. 

Together with the simple daisy and bridge classes introduced in Fig. [21 one can 
contemplate more complicated corrections. Bridges can be chained together in arbitrary 
combinations, as shown in Fig. IBl( a). or the component lines within a given bridge can 
themselves be joined together by other particles, as in Fig. IBlf b). Alternatively, bridges 
can be nested within each other to create rainbows — see Fig. IBl( c). In principle, a 
hierarchy of resummations (somewhat similar to the Balitsky hierarchy in QCD) is 
necessary to accommodate all these types of activity. 

Appendix B.l. Daisy diagrams 

Let us first consider the effect of daisies which dress bare Standard Model vertices. 
These are always momentum-independent and merely constitute a renormalization of 
whichever coupling constant sets the strength of the interaction at the vertex. For this 
reason they are relatively straightforward to deal with, and in the simplest situation we 
shall be able to resum their effect to all orders. If the daisies vary between different 
species of fermion, then the result would be an apparent species-dependent Fermi 
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Figure Bl. More complicated bridge-class diagrams. In addition to the simple bridge 
shown in Fig. HJc), one can use high-order vertices between two gauge bosons and an 
arbitrary number of dark energy quanta to "chain" any number of bridges together, as 
shown in (a). In (b), the component lines of a particular bridge are themselves joined 
together by virtual quanta of other species. (These could include dark energy particles, 
because the loops formed by such "joined bridges" would not already be included in 
the effective dark energy sector.) In (c), bridge diagrams are nested within each other 
to form so-called rainbow diagrams. 



constant, Gp- To prevent this occurring the fermion couphng function must be universal, 
which will be the case for conformal couplings. In what follows, we assume this to be 
the case. 

Consider Eq. ([2]) and expand the coupling functions B{x) and Bh{x) according to 

oo ^ 

^-^ nl 

n=0 

oo 

BniPx) ^ Yl -,BhM{^xT. (B.2) 

n=0 

where x = X + given that x is the expectation value of the dark energy scalar in the 
vacuum, and Bn {BH,n) are the Taylor coefficients of B (Bh) evaluated in this vacuum. 
We assume that Bq = Bhq = 1 and introduce a quantity 7„, defined by 

In terms of i?„ and 7„, the nth order interaction vertex takes the form 

f d^ki d^k2 dVi dVn ^4^/, , / , n 
^" = y (2^0^ (2^0^ (2^0^ ■ ■ ■ (2^(2vr) + k, + ^p,) 

-^W+{k,)W,-{k2)Sx{pi) ■ ■ ■ SxMv^'ih ■ k, - 7„M^) - k\k^]. (B.4) 

As an example, we will compute the simplest class of daisies which contribute to the 
interior of the propagator. The calculation of daisies which dress vertices with other 
species of fermion — or for the other gauge bosons — proceeds analogously. One finds that 
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the vacuum polarization with momentum transfer q, which arises from the daisy with n 
petals, can be written 



2n 



n 



V. 



2ti^k^ dK 1 



(2vr) 



Ml 



(B.5) 



{2n)\ 

where [n — 1)!! = {n — l)(n — 3) ■ ■ ■ 1 is the so-called "double factorial." In the special 
case of an exponential coupling, where = 1 for all n, and assuming that 7„ can be 
replaced by a constant 7, then it follows that all orders of daisies can be resummed to 
give 



1/4,3/4 
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A2 M2 
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A2 M2; 



3/4,5/4,3/2 



21^4 
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M2 



A2 



1 2 
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A2 M2 

A, 



(B.6) 



^ and < A 



where F{a; b\z) is the generalized hypergeometric function. For A < /3 
this resummation is dominated by its one-loop term. There will be extra terms in 
addition to Eq. flB.6l) which arise from interference between daisy diagrams and bridge 
diagrams (to be discussed in the next section). Although these may change the details 
of some numerical coefficients, they will not alter the momentum-independent character 
of the corrections. 

Eqs. (IB.5p - (lB.6p exhibit the general features which will recur in all diagrams we 
consider in this Appendix. A diagram with n dark energy lines can contribute at leading 
order in powers of (3 A, with this contribution coming from a region of phase space where 
all dark energy lines are carrying momenta of order A. This would seem to suggest that 
diagrams containing an arbitrary number of lines need to be accounted for in order to 
make reliable predictions. However, each line is also accompanied by a phase-space factor 
of 1/I67r2 (plus the combinatorical factor [n — l)!!/(2n)!) which leads to suppression of 
high-loop terms, so that counting powers of /?A alone does not give a proper accounting 
of the relative magnitude of adjacent terms in the loop expansion. 



Appendix B.2. Bridge diagrams 

Bridge diagrams are more complicated to handle. We proceed in two steps, first arguing 
that a similar phase-space suppression means that only the one-loop bridge need be 
considered, and not multi-loop or rainbow bridges. In a second step, we argue that 
chains of bridges can be ignored because they make a contribution at leading order in 
powers of /9A which is precisely momentum-independent. A contribution of this type can 
be absorbed into coupling constants and becomes unobservable. The one-loop bridge 
is itself momentum-independent at leading order in powers of /3A, so that the only 
contributions at this order which are not suppressed by a phase space factor of order 
~ 100 are the oblique corrections considered in ^ It will transpire that we expect 
corrections to the purely oblique analysis of §^|3H5]to occur at a relative order of roughly 
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l/Svr^ ^ 0.013 or better. We believe this is an acceptable precision at which to predict 
what can be observed at present and future particle colliders. It is again important that 
there is a universal fermion coupling function, in order that the effective Fermi constant 
Gp does not become species- dependent. 

First consider a multi-loop contribution to the vacuum polarization of any SU (2) 
gauge boson. At momentum transfer q, a calculation similar to those presented in ^ 
establishes that this can be represented in the form 
^„ ^Bl(3'\ f dH dVi dV„_i 



- Bin\ ^ ' J (27r)4 (27r)4 (27r)4 

P + M^-ierf + M^-ie rl_^ + - ie + - ie' ^ ' ' 

where R = q + i + YTjZi "^j-, ^ is the momentum carried by the exchanged gauge boson, 
and Pab{(^,q) satisfies 

Pab = Vab{i ■ q - InM^f - {iaqb + hqa){^ " q - InM^) + lJb{q^ + llMlr). (B.8) 

Eq. (1B.7P is suppressed by 2n powers of the coupling /3, but each scalar integral can 
contribute a power of A^. Since Pab ~ the C. integration can contribute terms of 
order A^. This would seem to imply that terms of order h?{(3Ky'^ could be present in 
the answer, but the correct conclusion depends on the relative magnitude of R^. Unlike 
the daisy diagrams or one-loop bridges, there can be some regions of phase space where 
£ ~ r j ~ A and |£ + Yl^=i ^il ~ 0; so that R ^ q. However, if the region of phase space 
in which this finely-tuned cancellation occurs shrinks with increasing cutoff faster than 
A^ then we can estimate the leading contribution by setting ~ A in Eq. (IB. 711 . In 
practice, the enhanced region of phase space is negligibly small. 

The n-loop bridge. To estimate the contribution of the n-loop bridge, we set R ~ nA 
and replace each factor such as d^rj/(27r)'^ by A^/Stt^. This choice for R is tantamount 
to assuming that the rj and i are randomly oriented, so that their cross terms 
approximately average to zero. This is likely to be a good approximation for large 
n but may fail for n ~ 0(1), so we will demonstrate explicitly that this procedure gives 
the correct answer for the 2-loop bridge. We consider this to be reasonable evidence 
that our estimate is reliable for all n. Proceeding in this way, it follows that the n-loop 
bridge makes a contribution to the vacuum polarization which is roughly equal to 

For the special case of an exponential coupling — for which Bn = 1 for all n — this can 
be resummed to give 



Kb^7Vabi2q' + ^'^M'^)^tL±_[ ] . (B.9) 



1 r d'^A'^ B'^A^ 

I^ab ^ -Vabi2q^ + 7'M^) <^ -^E + Ei( V^) - In ^ \ , (B.IO) 



4—--' ■ ' - -yy' ^ - ■ -V ^ 8^2 

where je ~ 0.577 is the Euler-Mascheroni constant and 'Ei{z) is the exponential integral. 
As in the case of the daisy diagrams this is dominated by its one-loop contribution, with 
the contribution of higher loops being suppressed by the phase-space factor l/Svr^. 
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Two-loop bridge. In the special case of the two-loop bridge, Eq. (1B.7P can be evaluated 
by the method of Feynman parameters, as in §U with the result that 

nl? = -f- [^j I dxdyd^ 5(1- ^-y-^) 

where F is defined by 

F = (1 - xy + Wf + XYq^ + xM^ + (1 - x)MJ, (B.ll) 
and the three quantities PF, X and Y satisfy 

W . (B.12) 

1 — X 



+ (B.14) 

Performing the I and r integrals and keeping only the leading term in powers of /3A, 
this vacuum polarization can be simplified to read 

ni? - -^^(2g^ + llMl)^^ [t^ t ^^^^^^ 5(1 - X - y - 



1 

X 



+ (W^ + 1 - x) In ■ 



X 



(B.15) 



Vr + 1 - x_ 

where by "~" we mean that this relationship is true up to terms of order /3^(/5A)^ 
which we have neglected. The possibility of enhanced regions of phase space where 
the loop momenta approximately cancel to leave an anomalously small propagator 
~ (rather than ~ has been replaced by the potential for large contributions 

from the Feynman parameter integrals. Indeed, inspection of Eq. (IB. Ill) shows when 
considering only the leading term in powers of A we might find a divergence roughly 
like J dx/ (1 — x)^. This would be finite when terms of all orders in A were included, but 
would manifest as an apparent divergence in the truncated series. If such a divergence 
appears, it should be regulated at a scale roughly of order M^/A^ where other terms in 
the perturbation theory become important, allowing enhanced phase space regions to 
appear. However, when these integrals are treated sufficiently carefully we find that 
no divergences occur and therefore that no enhanced regions of phase space exist. 
The integral can be evaluated by an adaptive Monte Carlo technique, and we find 
its numerical value to be roughly ^ 0.025. We conclude that Eq. (IB. 151) is as small in 
magnitude, or slightly smaller, than our estimate Eq. fIB.lOp . 



Chains of bridges. Now consider chains of bridges. We will first give an argument 
that the leading term in powers of /3A is entirely momentum independent for a one-loop 
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(a) k 



(b) 



Figure B2. Bridge diagrams. In (a), 2 — > 2' fermion scattering is dressed by a single 
dark energy bridge, in which a single dark energy particle is emitted or absorbed at 
the scattering vertices. In (b), a chain of bridges is exchanged. In principle, more 
complicated configurations exist in which the bridge is itself built out of rainbows of 
sub-bridges. We neglect these, since we anticipate that they will be suppressed by 
extra powers of the phase-space factor I/Stt^. 



bridge, before generalizing this to chains of arbitrary length. In each case, we use the 
results of the previous section to drop terms containing more than a single loop. 



The single bridge is shown in Fig. IB2f a). We assume that it describes a bridged 

fermion scattering process, between two fermion species which couple conformally to 
dark energy via the coupling function F{(3fx) and coupling f3f. (The generalization to 
different couplings and coupling functions is obvious.) It is easy to see that the effect of 
the bridge is the same as inclusion of an extra term in the propagator, corresponding to 



Clearly, the leading term in powers of /3A is momentum independent. 

The case of multiple bridges chained together is shown in Fig. IB2( b). One may 
wonder why it is necessary to consider such chains, since we have already argued that 
diagrams with a large number of loops are suppressed by powers of the phase-space 
factor Svr^. The reason is that to study corrections to the propagator, which potentially 
shift the location of its pole, one must resum enough diagrams to capture shifts in 
the physical mass of the particle. Such shifts are not captured at any finite order in 
perturbation theory. 

Exchange of the chained bridge with n internal vertices is equivalent to introducing 
an extra term in the propagator, which takes the form 





where A satisfies 




(B.17) 




ie (g + kf + M2 - ie 
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In the first contraction, the cubic term ~ cancels out. In the second contraction, 
the cubic term ~ cancels out. The same sequence recurs throughout the chain. This 
implies that the only way to bring each momentum integral to order ~ is consider only 
the quadratic term from each propagator, which multiplies the momentum-independent 
quantity 7iyM^. Accordingly, one can conclude that the resummed contribution at 
leading order in (3K behaves like an extra term in the propagator of the form 

^f?(/J/A)^,.J^-^. (B,18) 

Once again, this term is momentum-independent and contributes only to an 
unobservable shift — the same for all species — in the relevant masses and coupling 
constants. 
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